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Abstract. We consider two two-level atoms fixed at different positions, driven 
by a monochromatic laser field, and interacting collectively with the vacuum 
electromagnetic field. A Born-Markov-secular master equation is used to describe 
the dynamics of the two atoms and their steady-state is obtained analytically for 
two configurations, one in which the atoms are in equivalent positions and another in 
which they are not. The steady-state populations of the energy levels of the free atoms, 
entanglement, quantum discord and degree of mixed- ness are calculated analytically 
as a function of the laser field intensity and the distance between the two atoms. It 
is found that driving both atoms with the laser field is inefficient for the generation 
of steady-state correlations when they are in equivalent positions. On the contrary, 
inequivalent positions lead to the possibility of considerable steady-state entanglement 
and left/right quantum discord. It is shown that an X-state can be obtained for 
high laser field intensities for both configurations. This allows the comparison of two 
measures of quantum discord. The behaviour and relationships between correlations 
are studied and several limiting cases are investigated. 
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1. Introduction 

For the past two decades there has been an ever growing interest in quantifying and 
characterizing the correlations in the states of quantum systems. This was initially 
provoked by the realization that quantum correlations known as entanglement could be 
harnessed as a resource to bring significant advantage for computing and information 
processing [1]. Now, one of the main interests is to identify which correlations are 
responsible for these advantages. Moreover, it has been recognized that quantifying 
correlations in quantum systems is a difficult task and, in general, several measures are 
needed to capture all of their subtleties. 

In order to understand correlations, one generally starts with bipartite systems. 
Consider two quantum systems A and B. The state of the composite system A + B 
can be described by a density operator p^g which contains both classical and quantum 
correlations. A widely accepted measure of the total correlations in p_4s is the quantum 
mutual information. It is defined as 

Hpab) = S{pj^) + S{pb) - S{pab) , (1) 

where p_A (pb) is the density operator of A (B) and S{p) = — Tr[plog2(p)] is the von 
Neumann entropy of p (the trace is taken in the state space where p is a density 
operator). It has been found that I{pab) is the distance (as measured by the relative 
entropy) of pj^q to its closest product state p_4 ® pg [2J. It has also been shown that 
I{Pab) measures the asymptotically minimal amount of local noise one has to add to 
turn pab into a product state [3]. Finally, I{pab) is the maximum amount of information 
that Alice can send secretly to Bob using p^B as a one-time pad [4]. 

Once the total correlations in p^B have been quantified, it is natural to ask whether 
these can be clearly divided into a classical and quantum part. Several axioms have 
been proposed as requirements of a measure of the classical correlations C'^\pab) in 
Pab These consist in being zero for product states, being invariant under local 
unitary transformations and non-increasing under local operations, and being equal to 
S{,Pa) = S^Pb) for pure states pab- Also, it has been pointed out that a measure of 
the classical correlations should quantify the correlation between A and B instead of a 
property of only one of them, that is, it should be symmetric under interchange of A and 
B. Taking these axioms as a basis, the following measures for the classical correlations 
in Pab have been proposed [Sj: 

Cb{pab) = S{pa) -T^in{B,} ^Pi5(p^) , 

i 

CUpab) = S{pb) -min{^,} Y^P'^^W , (2) 

i 

where the minimum is taken over all sets {-Bj} ({Aj}) of one-dimensional orthogonal 
projectors that sum up to the identity and that constitute measurements performed 
only on B {A). Also, 

p\ = -i:iB\{l®Bi)pAB{l®B\)\ , (3) 
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with 



Pi = Tr^B 



{I®B,)pM^®bI) , (4) 



is the density operator of A after obtaining the result associated with Bi in a 
measurement of B (similar equations hold for pg and p'j). These measures satisfy 
the axioms mentioned above and are equal to zero if and only if p^g = pa ® Pb- 
Hence, it is in accordance with the accepted idea that only product states are devoid of 
correlations. Nevertheless, these measures are dependent on which system is measured, 
that is, C^(p^b) 7^ C^^p^ab) in general. Thus it depends on the properties of each 
subsystem. 

More recently, other measures and notions of classical correlations have been 
proposed [2], [3]. In particular, it has been proposed that the classical correlations 
should be quantified by the minimum distance (as measured by the relative entropy) 
between a classically correlated state associated with p_4B and a product state [2] . This 
last proposal has the advantage of placing all correlations (total, quantum, and classical) 
on an equal footing that allows a direct comparison between all of them [2]. 

Given that the classical part of the correlations is not yet clearly quantified, let 
us now turn to quantum correlations. Quantum states p_ab are normally divided into 
separable or entangled. Let us remember that p^g is separable if it can be expressed in 
the form 



Pab 



^PjpA,j® PB,j ^ (5) 



with p_A,j {pbj) density operators of A {B), and pj G [0, 1] such that J2jPj — 1- Pab 
is not a separable state, then it is an entangled one. It was thought that entanglement 
embodied all the quantum correlations in p^g, and that separable states were purely 
classical. Nevertheless, it has been realized that entanglement is not the only aspect 
of quantum correlations, since some separable states may still present non-classical 
correlations [2]-[H]- One way to measure the non-classicality of the correlations in p_4g 
is to use the quantum discord [6]. This quantity is defined to be the difference between 
the quantum mutual information ([T]) and the correlations ([2]): 

D%Ab) = IiPAB)-Ci{PAB) , 

D^iPAB) = I{pab)-CI{pab) , (6) 

Dj^{Pab) {D^{pab)) is usually referred to as the left (right) quantum discord. It has 
been shown that D'^{pab), D^{pab) are always non-negative [6]. In fact, -D^(p^g) = 
if and only if 

P = ® Pfc , (7) 

k 

where {I'l/'fc)} is an orthonormal basis of A, pk are density operators of B, and pk are 
non- negative numbers such that J2kPk = 1 [6j. Similarly, Dq{pab) = if and only if 

p = y^^PkPk ® \i^k){i^k\ , (8) 



Steady- state quantum correlations 



4 



where {iV'fc)} is an orthonormal basis of B, pk are density operators of A, and pk are 
non-negative numbers such that YlkPk — 1- 

In general -D5(p^s) 7^ -Dg (p^e)- The quantum discord D^^pj^s) can be interpreted 
to be a measure of the information of A in p_4g that cannot be extracted by performing 
measurements only on ,B [6]. Hence, if -D^(p^b) is large, a lot of information is destroyed 
by any measurement on B alone; if D^^py^^ig) is small, almost all the information about 
A in p_4g can be obtained by measurements only on |6]. 

Once quantum discord is recognized as a measure of non-classical correlations, it is 
natural to ask if it measures the same correlations as entanglement. Given that some 
separable mixed states have non-zero quantum discord and separable states by definition 
do not have entanglement, it is concluded that entanglement and quantum discord are in 
general different quantities. Nevertheless, they do coincide when pyig is a pure state or 
a mixture of Bell states [7]. Furthermore, it has been found that in general the quantum 
discord of p^g is not simply the sum of some measure of the entanglement in p_4g and 
some other non-classical correlation [7], [8]. 

The study of the dynamics of quantum correlations is intimately linked to the 
study of open quantum systems, that is, quantum systems which are coupled to other 
quantum systems called environment [9|. When the environment has many degrees of 
freedom (a reservoir), it normally induces a loss of coherence in the system of interest 
that ultimately has a detrimental effect on the correlations present in it. For example, 
entanglement may decay exponentially with time or may vanish completely in a finite 
time (called entanglement sudden death) [10]- [12]. Moreover, the environment can also 
affect the degree of mixed-ness of the state of the system of interest [9], and this in 
turn also affects the entanglement since it is known that a mixed state cannot have an 
arbitrary degree of entanglement [13] . 

Recently, the open system dynamics of both the quantum discord and classical 
correlations have started to be studied [ll]-[I7]. In particular it has been found that 
under Markovian environments quantum discord may be more robust than entanglement 
[T5] . For non-Markovian environments and independent reservoirs for each part of a 
bipartite system, the quantum discord may vanish only at discrete instants whereas 
the entanglement can disappear during a finite time interval. For a common reservoir, 
quantum discord and entanglement can behave very differently with sudden birth of the 
former but not of the latter [T7] . 

It is the purpose of this article to study the entanglement, quantum discord, classical 
correlations, and degree of mixed-ness in the steady-state of the following open quantum 
system: two two- level atoms (qubits in the jargon of quantum information) at fixed 
positions driven by a monochromatic laser field and interacting collectively with all the 
modes of the quantum electromagnetic field. The latter will be assumed to be in the 
vacuum state. In the notation used above, A will be one of the atoms (say, the atom at 
position ri) and B will be the other atom (say, the atom at position r2). Here we also 
have a third party C playing a decisive role in the dynamics: the vacuum electromagnetic 
field which we will consider as a reservoir. The system is open because we are interested 
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in studying only the subsystem A + B of the complete system A + B + C. Since the 
two atoms are interacting collectively with the reservoir, C acts as a medium that can 
allow quantum correlations to be formed between the two atoms for some time. In order 
that these correlations have the possibility of being long-lived and not being ultimately 
destroyed by the reservoir, the two atoms will be driven by a laser field. 

We will investigate in detail two configurations for the two atoms. In the first, one 
of the atoms will be placed in a position where the laser electric field is zero, while the 
other atom will located at a position where it is not zero. In the second, each atom 
will be placed in a position where the driving electric field takes the same value. These 
two configurations will allow us to study how the correlations change when the atoms 
are in equivalent /non-equivalent positions and to determine what is the effect when 
both atoms are being driven by the laser field. It is important to mention that we will 
neglect the coherent dipole-dipole interaction between the two atoms that results of the 
collective interaction with the reservoir. We will only keep the dissipative collective 
interaction with the reservoir. This will allow us to calculate a steady-state density 
operator of the two atoms analytically, as well as, all the aforementioned correlations 
in it. Furthermore, this will allow us to identify which effects are due solely to the 
dissipative interaction. In future work we will add the dipole-dipole interaction. 

Using the same system as ours, but without the driving field, the dynamical 
generation of entanglement between two atoms due to the collective interaction with 
the reservoir has been recently studied [I8],[19]. It was found that, as a result of the 
interaction through the reservoir, the system develops non-negligible entanglement (as 
measured by the concurrence) for a period of time if the atoms are close enough and if 
initially one is in the excited state and the other is in the ground state. Also, it was 
reported that entanglement initially present in the system of two atoms is more robust 
when the two atoms are close when compared to the case where the atoms are far apart. 
The entanglement present in a system of two driven non-identical two-level atoms in 
a special configuration (namely one atom is located at a node while the other atom is 
placed at an antinode of a driving laser field with a standing wave cosine structure) has 
also been studied [20]. The steady-state density operator of the system was obtained 
numerically and the concurrence was evaluated. It was found that the two atom system 
decays to a stationary entangled state only when the Rabi frequency equals the difference 
between the two atoms' transition frequencies. 

The present article is organized as follows. In Section II we summarize some 
results on measures of correlations for two qubit systems. In Section III the system 
of interest is described and the master equation governing the dynamics of the two 
atoms is established. In Section IV, the steady-state density operator is calculated 
analytically for the case of only one atom being driven by the laser field. The populations 
of the eigenstates of the free Hamiltonian of the two atoms, the degree of entanglement, 
the quantum discord, and the degree of mixed-ness of the two atoms are evaluated as 
functions of the laser field intensity and of the distance between the atoms. In Section 
V, the steady-state density operator of the two atoms is calculated analytically for the 
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case of equal coupling strengths. The same analysis of Section IV is done for this 
new configuration. Furthermore, the results of both configurations are compared. The 
conclusions are given in Section VI. 



2. Quantifying correlations in two qubit systems 



In this article we will be considering two two-level atoms (qubits) which we will number 
by 1 and 2. In the following sections these labels will correspond to the atom at 
position ri and to the atom at position r2, respectively. In terms of the notation of 
the Introduction, 1 will replace A, while 2 will replace B. We will now show how to 
calculate all correlations by simple formulas. We will only give the algorithms and refer 
the interested reader to the original articles for the proofs. 

The kets \j : +) and \j : — ) will denote the excited and ground states of the jth 
atom (j = 1,2), respectively. In the following we will be making constant use of the 
triplet-singlet basis 

B = { |1,0), |0,0) } , (9) 

for the state space of the two atoms: 



|1,1) = |1 : +) ® |2 : +) , 
|l,0) = i=( |1:+)®|2: 
|1,-1) = |1:-)®|2:-) 
|0,0) = i=(|l:+)®|2: 



-) + |1 : -) ® |2 : +) ) 



-)-|l:-)®|2:+) ) 



We shall also use the usual tensor product basis 
where 



B' = { 


+,+), 1+,- 




1+,+) = 


|1 : +) ® |2 


+) 


1+,-) = 


|1 : +) ® |2 


-) 


1- +) = 


|1 : -) ® |2 


+) 




|1 : -) ® |2 


:-) 



-) } 



(10) 



(12) 

Furthermore, we will denote the density operator of the two atoms by pi2. 

We are interested in quantifying the degree of entanglement of the system of two 
atoms. We will use the concurrence C, which can be calculated as |2T] : 

(13) 



C = max{0, V Ai — V A2 — V A3 — V A4} , 
where Ai > A2 > A3 > A4 are the eigenvalues of the matrix 

M= [pi2]B'{(ry(g)ay)[pi2]*B'i(^y'S)cry) . (14) 

Here ay is the well-known Pauli matrix, [pi2]b' is the matrix representation of pu in 
the basis B', and [pi2]b' is the element-wise complex conjugate of the density matrix 
[Pi2]b'- Note that any basis of the state space of the two atoms related to the tensor 
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product basis by an orthogonal change of coordinates matrix can be used for the matrix 
representations (for example, B could be used instead of B'). The concurrence takes 
values between and 1. It is 1 when the atoms are in a maximally entangled state, 
while it is zero when the atoms are in a separable state. 

Evaluation of the quantum discord given by (EI) in general requires considerable 
numerical minimization. Although a method to calculate easily the classical correlation 
and quantum discord for a general two-qubit X-state has been proposed [T] , it has been 
found to be unreliable because it does not take into account all of the restrictions in the 
minimization problem [22]. Hence, we will show here how to calculate numerically the 
right quantum discord. 

The first step is to express pi2 as a linear combination of Pauli operators and tensor 
products of Pauli operators: 

Pi2 = ^ ^II + x-d-i + y-d-2+ ^TijCTiiOaaj ^ . (15) 

Here I is the identity operator, aj is the vector of Pauli operators of atom j and an ® o"2j 
is the tensor product of the Pauli operator an of atom 1 and the Pauli operator a2j of 
atom 2. Notice that x is the Bloch vector of the density operator pi = Tr2[pi2] of atom 
1, while y is the Bloch vector of the density operator p2 = Tri[pi2] of atom 2. 

The next step is to see what is the structure of a one-dimensional orthogonal 
projector Pq in the state space of atom 2: 



a 



P' 1b, = I J,.. 1 _ „ j . (16) 

where B2 = {|2 : +), |2 : — )}, [Po]b2 is the matrix representation of Pq with respect to 
the basis B2 of the state space of atom 2, a e [0, 1], (f) E [0, 2tt) and b = a/ a{l — a). 
With (|T6|) at hand, we notice that 

Poa^rPo = 26 ■ cos(0)Po , 
Poa,Po = 26-sin(0)Po, 

Poa^Po = (2a-l)-Po . (17) 
We now define 

7 = (71, 72, 73) = (26-cos(0), 26-sin(0), 2a - 1) . (18) 

Using f|T5|) and f|T7|) we now see how pi2 transforms when a measurement on atom 
2 is performed and the result associated with Pq is obtained: 

-(I ® Po)pi2(I ® Po) = f pl + V iy,a, I ® Po , (19) 

Po 4po V tr y 

= PPo ® Po , (20) 



where 



Po = Tr[(I®Po)pi2(I®Po)] = ^p. (21) 
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Here we have introduced the quantities 



/i 



y -7, 



+ ^ = (^^1^1^2,1^3) ■ (22) 

i=i 

Note that pp^ is a density operator of atom 1. Since Pq is an orthogonal projector in 
the state space of atom 2 it follows that 

(23) 



pPo = Tr2<^ -(I®Po)pi2(I®Po) 
We now calculate the von Neumann entropy of pp^ 



S ( PPo) 



— ) log2 
p 



1 f I'^l 



1 - 



r I 
P 



(24) 



With results fll5p -f l24p at hand we are all set to determine a formula for the right 
classical correlations Cf[pi2)- Let {!'?/'), |0)} be a complete set of one-dimensional 
orthogonal projectors (that is, an orthogonal basis) for the state space of atom 2. 
Using the basis B2 = {|2 : +), |2 : — )} it is straightfoward to show that the matrix 
representations of and |0) have the form 

' |2 



I ape** 



1 — a 



1 - 


a 


2 






a|2 



(25) 



<l> in all 



where |a| G [0,1], $ G [0, 27r). Notice that if we take a = |ap and (j) 
the results for Pq we will get the corresponding results for Similarly, taking 

a = 1 — |ap and = $ + tt we get the corresponding results for If we define 

Pip = PPo when we identify Pq with \ilj){ip\, and we take = pp^ when we identify Pq 
with 10) (01, we find using ([T5])-([2l]) that 



P^S{p^) + P4,S{p^) , 

loga 
log2 
log2 
log2 



p+ 


(- 


1^^+ 




4 




/i+ 


P+ 


(- 


\^+ 


4 




P+ 


P- 


(- 


z/_ 


4 






P- 




z/_ 


4 




P- 



p+ 



+ 1 



1 + 
1 - 



p+ 

|z/_ I 

|z/_ I 
P- 



(26) 



where we have defined 

li> = (7vi' 7^2, 7V3) 
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7^1 = 2 ■ cos(<l>)|a|A/l - , 
7^2 = 2 ■ sen($)|a|v^r^^ , 
7^3 = 2|«|^ - 1 , 

u±i = Xi ± T(i) ■ 7^ , 

/x± = 1 ± y ■ 7V, . (27) 

Also, T(j) is the i-th row of the matrix T whose component i,j is Tij. 

To calculate the classical correlation Cf{pi2) we have to minimize (!26|) over all 
complete sets of one-dimensional orthogonal projectors { |0) }. From ( 125|1 we note 
that the set { l?/^), |0) } is determined by two real parameters: |a| G [0, 1], $ G [0, 27r). 
Then Cf{pi2) can be calculated as follows: 

C2\pi2) = S{pi) - min|c,|e[o,i], $e[o,2^] [Pi^S{p^^) + p^S{p^)] , (28) 

where PiijS{ptp) + p<f,S{p^) is defined in fl26|) . We note that the minimum always 
exists because the minimization process is done over a compact set. Also notice that 
$ G [0, 27r) was extended to $ G [0, 27r]. This was done to have a compact set and is of 
no consequence since we are interested in the absolute minimum of a periodic function of 
period 27r. Also notice that once C2{pi2) has been calculated, the right quantum discord 
D^(pi2) follows easily by substracting Cf{pi2) from the quantum mutual information 

I{Pl2). 

Now suppose that we make measurements only on the atom at position ri, that is, 
we want to calculate the left classical correlations (pi2). Instead of going through the 
whole process of calculating conditions analogous to those in fl26l) - fl28l) . we only need to 
find the matrix representation of the density operator pi2 with respect to the basis B" 
in which the order of atoms 1 and 2 has been interchanged 

B" = I |2 : +) ® |1 : +), |2 : +) ® |1 : -), 

|2 : -) ® |1 : +), |2 : -) ® |1 : -) } , (29) 

and apply the procedure described above in fl2B]) - fl25]) to obtain the quantum discord 

Di{pi2) and the classical correlations Cf{pi2). It is easy to find such matrix 

representation. Suppose that the density operator pi2 has the following matrix 
representation with respect to B' 

/ Rii Ri2 Ri3 -Ri4 \ 

i?2l R22 R23 R24 
R31 R32 R33 R3A 
Y i?4l R42 R43 -R44 / 



[p 



12 B' 



(30) 
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Then the matrix representation of pi2 with respect to B" has the following form: 



[P 



12 B" 



( Ru 


Rl3 


R12 


R14, 


\ 


R31 


R33 


R32 


-R34 




R21 


R23 


R22 


-R24 




\ -R41 


-R43 


Ra2 


-R44 


/ 



(31) 



Note that [pi2]b" is obtained by interchanging columns 2 and 3 and rows 2 and 3 of 

[P12]b'- 

For a general two-qubit density matrix the quantum discord has to be calculated 
numerically by a procedure such as the one described above. Therefore, several 
alternative measures have been proposed. One in particular calculates the distance 
of P12 to the set of zero discord states Qq given by ([7]) if measurements are made on 
atom 1 or given by (IHl) if measurements are carried out on atom 2 ^23j. It was found 
that j23] : 

1 



12 



mm 



ollPi2 - Xl 



|2 

If 



(xx^ + 



|T|||- A;- 



max 



) 



(32) 



Here x is a real column vector whose three components are given by Xj = Tr(pi20"ii 01), 
T is a 3 X 3 real matrix whose components are given by Tjj = Tr(pi2crij ® o'2j), and 
fcniax is the largest eigenvalue of the matrix K = xx-^ + TT^. Note that aij is the j-th 
Pauli matrix of atom i. Also, the distance is measured using the usual Hilbert-Schmidt- 
Frobenius norm: 



A\\l = Ti{A^A) , 



(33) 



with A a linear operator in the state space of the system (it could also be a square 
matrix). D{ (pu) is called the geometric measure of left discord. Using (l30l) and (l3Til 
one can determine the geometric measure of right quantum discord D2 {pu)- 

To quantify the degree of mixed-ness we will use the linear entropy Sl (9] defined 

as 



Sl{P12] 

Recall that Sl{pi2) = 
mixed state. 



= 1-Tr[py • 

if P12 is a pure state, while Sl{pi2) 



(34) 

3/4 if pi2 is a maximum 



3. Two driven qubits collectively interacting with a vacuum reservoir 

We consider two identical two-level atoms with transition frequency ua_ at fixed positions 
ri and r2, driven by a classical monochromatic laser field, and interacting with all the 
modes of the quantum electromagnetic field. In the following we will refer to the latter 
as the reservoir and, in some occasions, to the atom at position as atom j. 

The Hamiltonian of the system in the long-wavelength approximation and in the 
electric dipole representation is 

H = Ho + V + VAiit) (35) 
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where Ho is the free hamiltonian of the two atoms and of the reservoir 

Ho = ^CT3 + ^hwj (a]aj + , (36) 

V is the electric dipole interaction between the two atoms and the electric field E(r) of 
the reservoir 

V = -di- E(ri) - d2 • E(r2) , (37) 

and VAhit) is the electric dipole interaction in the rotating- wave-approximation between 
the two atoms and the classical monochromatic electric field Ei(r,t) of frequency ul 

2 

VALit) = -5^nG(r,)(a+,e-^-^* + a_,e-^*) . (38) 
i=i 

Recall that \j : +) and \j : — ) are the excited and ground states of the jth atom 
(j = 1,2), respectively, and 0-3 = (0-31 + (T32) with 0-3^ = \j : : +| - \j : : -| 

is the inversion operator. Furthermore, dj is the dipole operator of the jth atom and, 
since the atoms are identical and have two levels, can be expressed as 

dj = doicr+j + d*o^a_j , (39) 

where a±j = \j : ±)(j : =f| are the transition operators for the j-th atom and 
doi = (1 : +|di|l : — ). Notice that the dipoles have been taken to be parallel to 
each other. The quantum electric field E(r) at position r is given by its expansion in 
terms of plane waves 



where V is the quantization volume, aj(ap is the annihilation (creation) operator of a 
photon in mode j, and is a sum over all the modes of the quantum electromagnetic 
field. The wave and polarization vectors of mode j are kj and ij, respectively, while 
Uj = ckj is the corresponding angular frequency. 

We have taken the driving electric field Ei(r, t) of the form 

Ei(r,t) = g{r) {SLe-'^-' + Sle'^-') (41) 

with g{r) a real- valued function describing the spatial structure of the field and Sl 
a constant complex vector which contains the polarization of the electric field. Also, 
G{r) = 5f(r)(doi • Sl)/^ is assumed to be a real quantity, and in the following we 
shall denote G{rj) by Gj. For example, f HTj) could describe a stationary wave with 
polarization independent of position by taking g{r) = cos(k/^ ■ r). 

We will denote the density operator of the complete system (two atoms plus 
quantum electromagnetic field) by p{t), while pi2(t) will denote the density operator 
of the two atoms. Recall that pi2(t) is the reduced density operator obtained by tracing 
p{t) over the reservoir degrees of freedom. Also, we assume that the initial state of the 
system is of the form 

p{t = 0) = pi2(t = 0)® |0)(0| , (42) 
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that is, the system will initially be in a product state in which the reservoir is in the 
vacuum state. 

When the interaction with the vacuum reservoir is neglected (V^ is zero in ( 135|) ). 
Puit) is determined by von Neumann's equation: 



-(X3 + I/alW, Pl2{t) 



(43) 



On the other hand, when there is no driving field (Val(^) is zero in (1351) ). the dynamics of 
the density operator pi2{t) of the two atoms can be described by a Born-Markov-Secular 
master equation [24j: 



d 
dt 



Pl2{t) 



I 



hjjj , 



+ nPl2{t)] ■ 



(44) 



Here Hls is the Lamb-shift Hamiltonian which can be expressed as a sum of two parts: 
one that describes an effective coherent dipole-dipole interaction hdd between the two 
atoms and another that represents a shift eis of the energy levels of the free Hamiltonian 
of the two atoms [21]. The dissipator V is given by 



V{p) = V,{p) + V,{p) + V;{p) , (45) 
where Vj is the dissipator for the jth atom interacting with the reservoir (j = 1, 2) 



(j 



1,2), 



(46) 



and ©3 is the part that describes the dissipative collective interaction of the two atoms 
with the reservoir 

3 
2' 

(47) 



) / 1 

j7iF(|ri - r2|)(^a_ipa+2 - - {a+2(7-i,p} 

1 . , \ 



+ O--2PO-+1 - - {0-+i0-_2, p} 



Here {■, ■} is the anti-commutator and 71 is equal to the spontaneous emission rate of a 
two-level atom interacting with all the modes of the electromagnetic field 

1 4|doipa;i 
47reo 3hc^ 
The function F(|ri — r2|) is defined by 

d^, 



F(|ri-r2|) = 



loil 



sm X 



X 



+ 3 



loil 



cos X — 



sm(x) 



(4^ 



with 



X 



— 1^1 



^2\ 



and d^ the square of the norm of the projection of doi onto the plane perpendicular 
to ri — r2. In the following we will denote F{\ri — r2|) by F12 [27]. We note that Fu 
is an oscillatory function whose absolute maximum and minimum values are 2/3 and 
—0.2237, and that they occur only at the points — t'2\/c = (independent of the 

value of d\) and w^jri — r2|/c = 4.233 (with d^ = |doip), respectively. 
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Applying the approximation of independent rates of variation the master 
equation for puit) which takes into account the interaction of the two atoms with 
the vacuum reservoir and with the classical field is given by {t > 0) 

hoJA 



d , . i 
dtP''^'^ = —h 



-as + HLs + VAL{t), pi2(t) 



V[p,2{t)] . (49) 



2 

In the following we will assume that the energy shift e^s has been incorporated in Huja 
and we will neglect the coherent dipole-dipole interaction h^d between the two atoms 
which arises from the collective interaction with the reservoir. In general, h^d is negligible 
for distant atoms and for particular distances between the atoms |20] , but it can become 
important when the atoms are very close [21]. Later on, we will make explicit when the 
effects of hdd can be important. Furthermore, these will be taken into account in future 
work. 

Passing to the interaction picture (IP) defined by the unitary transformation 
Uo(t,0) = exp[—iuAO'3t/2] we obtain the master equation that will be taken as the 
model under study for the rest of the article: 

j^Puit) = KlW,Pi2W] +'D[p,2{t)] (t > 0) , (50) 



where 



Vidt) = - Y.hG{r,){a^,e-^'-' + a.,e''-') , (51) 
i=i 

5l = ujl — ua, and pi2{t) is now the IP density operator of the two atoms. 

To study the dynamics of the system we will calculate the steady-state solution pf^ 
of fl50l) for two special configurations which will be analyzed in the following. Recall 
that pf2 is a steady-state solution of ( l50l) if pf^ is not explicitly time dependent and 

-'j^[vL{t\plI]+V[pl^]=^. (52) 

We note that in general f l50|) has no steady-state solutions if 5l 7^ 0, since density 
operators pf^ that satisfy ( 152|) have non-diagonal time dependent elements if 5l 7^ 0. 
This will be shown explicitly in the following. Also notice that other steady-state 
solutions may appear in other interaction pictures. 



4. G(r2) = 

In this part we will assume that one of the atoms is fixed in a position where the classical 
electric field (14T|) is zero, while the other atom is fixed in a position where it is not zero. 
Therefore, we take G2 = G{t2) = 0. Notice that increasing (decreasing) the intensity 
of the laser field increases (decreases) \Gi\ = |G'(ri)| (see the definition of Gi following 
f l4T|) ). Therefore, Gi can be made to vary by increasing or decreasing the intensity of 
the electric field ( 14T]) . Furthermore, the distance |ri — between the atoms can also 
be varied independently of Gi. For example, one of the atoms could be placed at a 
node of a stationary wave to have G2 = 0, the other atom could be placed anywhere 
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else (except at a position where Gi = 0), and the intensity of the field could be varied 
to have \Gi\ = |G'(ri)| take on any positive value. 

It is found that the populations of the solution pfj of ( l52l) in the triplet-singlet 
basis (fTOj) are given by 



(l,l|pf2^|l,l) 

(l,0|pfj|l,0) 

(1,-I|pf2'^|l,-1) 

(0,0|pf2^|0,0) 



^[l6F^,Gl + 9{4-F^,)F^,Gt] 
- [ fi^Gl + iy_Gt + V-Gl ] , 



1 - (1, l|pfj|l, 1) - (l,0|p?2^ |l,0) - (0,0|pt2' |0,0) 



ST\ 



- [ p+G? + u^Gl + r^^Gl ] 



(53) 



where 



Gi 



Gi/7i , 

G^(2048 + 1152F{^2) + ^^(2560 + 5184^^2 - 1296Fi^2) 
+ 01(864: + 144^1^2 + 486Fi^2) 
+ ^(2 + 3Fi2)^(2- 3^12)^(4-^^2) , 



32 * (16 ± I2F12 + 9F^2) 
18 * (32 ± 48Fi2 + 68^2 
9 * (2 ± 3Fi2)'(4 - F^^) 



(54) 



On the other hand, the coherences of pfj are given by 

.9^2 



(0,0|pt2'|l,l) =te 



K 



GiF,2 3(2 + 3Fi2)(4-Fj 



(l,l|pf2^|0,0) 

(l,0|pf2'^|l,l) 

(l,l|pf2^|l,0) 
(1,-I|pf2^|l,0) 

(l,0|pf2^|l,-l) 

(1,-I|pf2^|0,0) 



+ 8(4 + 2Fi2 - 3F^2)Gl 
(0,0|pf2^|l,l)* , 



^s.t^Q3F J 3(2- 3F,2){4-F^,] 



K 



8(4 - 2Fi2 - 3F^^)Gl 
(l,0|pf2^|l,l)% 

32(32 - I2F12 - 18^1^2 + '^7F^2)Gl 
+ 36(32 - 24Fi2 - AF^^ + 54F'^2 - ^Fi2)Gl 
+ 9(2-3Fi2)'(2 + 3Fi2)(4-Fj 



4/t 



(1,-I|pf2^|l,0)* 



4:K 



32(32 + 12F] 



12 



12) 



18F^, 



27Fl,)G{ 



+ 36(32 + 24Fi2 - 4F/2 - 54Fi^2 + 9^12)^1 
+ 9(2 + 3Fi2)'(2 - 3Fi2)(4 - F^^] 



(1,-I|pf2^|0,0)* 



(0,0|pf2^|l,-l) 
(1,-I|pf2^|l,l) 

(I,l|pf2^|l,-l) = (l,-l|pf2^|l,l)* 



2k 



e'^^^'Fi2Gl - 256Gt - 9(64)^^2^? 



9(16 - 40^1^2 + 9F] 



12; 
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(0,0|pt2 11,0) =— 512Gt + 3Q{lQ + lQF^'2-QFi2)Gi 

K L 

+ 9(16 - AOF^, + 9F^,) ] , 
(l,0|pfj|0,0) = (0,0|pfj|l,0)* . (55) 

Note that some coherences of pfj depend exphcitly on time. Hence, in the IP we are 
working, a steady-state solution (|52l) exists only for the resonant case (5l = 0). The 
steady-state solution in fl53l) and fl55|) would also be approximately valid for small enough 
detuning 5l, for instance for \S-i,t\ <^ 1. In the following we will restrict to the case where 
the (exact) steady-state solution exists, that is, we will take 5l = 0. 

Using the following argument it can be shown that the density operator puit) for 
the two atoms tends to the state pf^ above for any initial state Pi2(0) if 6*2 = and 
there is resonance (5l = 0). To that end, we take G2 = and 5l = 0, and we choose 
the triplet-singlet basis ffTU]) for the state space of the two atoms. (ISUl) with the initial 
condition pi2(0) can then be reexpressed as an initial value problem (IVP) of the form 

x(t) = Ax(t) , x(0) = xo , (56) 

where x(t) is a 16 component column complex vector associated with pi2{t), A is a 
16 X 16 constant complex matrix, and Xq is a 16 component column complex vector 
defined by pi2(0). The IVP in ( l56l) has a unique solution of the form [26] 

x(t) = cie-^^Vi(t) + ... + ci6e-^^«Vi6(t) , (57) 

where —Xj are the complex eigenvalues of A, Vj{t) are vectors whose components are 
polynomials in t of degree less than the algebraic multiplicity of the corresponding 
eigenvalue —Xj, and Cj are complex constants determined by xq. When an eigenvalue 
—Xj has an algebraic multiplicity equal to its geometric multiplicity, the corresponding 
vectors vj reduce to linearly independent eigenvectors associated with —Xj. It can be 
shown that zero is an eigenvalue with algebraic multiplicity equal to one. If necessary 
rearranging the order, we take Ai = 0. Hence, vi is an eigenvector of A associated with 
Ai = 0. Furthermore, the rest of the eigenvalues have negative real part (that is, Re(Aj) 
> for j 7^ 1). It then follows that 

x(t) cie~^iVi(t) = civi ift^+oo. (58) 

Now CiVi is associated with pf|" above. Hence, for any initial condition puit) tends to 
pf^ as t — > +00. 

There are several limiting cases of interest. First, if Gi — ?■ 0, then it is seen from 
fl53|) and fl55|) that pf^" — )■ |1, — 1) (1, — 1|. This result is expected since the atoms are 
located at different positions and, without the driving field, the reservoir ultimately 
leaves the two atoms in their respective ground states. 

Another limiting case of more interest occurs when the classical electric field is very 
intense (|Gi| — )► +00). From (I53l) - fl55l) we see that the matrix representation of pf^ in 
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the tensor product basis B' fll2p takes the form of an X-state 



STl 
12 Jb' 



1^12 







v 






3Fi: 



9p2 
4 12 



16+9F{. 
3Fi2 





9 p2 
4-'^12 

16+9^^2 



16+9F/. 











3Fi2 
T6+9F^ 




9 7712 

- — 

16+9R2. 



(59) 



In this case the corresponding populations of the triplet-singlet basis B ( ITOl) take the 
form 

9 



:i,l|pf2^|l, 



1) 



4-^12 



(l,0|pfj|l,0) 

(1,-I|pf2^|l,-1) 

(0,0|pfj|0,0) 



16 + 9F22 ' 



1 

4 

1 

2 
1 



16 + 9F2 



9 f2 
4 -'^12 



12 



16 + 9F^^ 



+ 



3Fi 



12 



16 + 9F^^ 



(60) 



Noting that F12 takes on its absolute maximum value of 2/3 when |ri — = 
(independent of the value of d'jj and its absolute minimum value of —0.2237 when 
|ri — r2| = 4.233 and d^ = |doip, we can estabhsh bounds on the populations in 
f l60|) . The populations of the states |1,1) and |0,0) cannot grow more than 1/20 and 
1/4+1/10, respectively, and these maximum values are achieved when the atoms are very 
close together and the electric field is very intense. At the same time the populations of 
the states |1, — 1) and |1, 0) take on their minimum values of 1/2 — 1/20 and 1/4 — 1/10. 
We note that the dissipative collective interaction of the two atoms with the vacuum 
reservoir is responsible for the non-zero population of the state |1,1) and that this 
interaction is not so effective, since this population can only grow to a relatively small 
value of 1/20 for high laser field intensities. 

Now lets take the limit of the two atoms very far apart (|ri — — > +00). From 
( 153|1 and f l55|) we get the following density matrix in the tensor product basis B' f lT2|) : 



(0 











\ 










„• 2Gi 




1+SG'i 


^ 1+8G^ 




















2Gi 





l+4Gf 


/ 


^ 1+8G^ 


l+8Gf 



(61) 



Note that in this case the system behaves as if the two atoms were interacting with 
independent reservoirs and, since only one atom is being driven by the laser field, the 
population of the state |1, 1) = |+, +) is zero. 

A final limiting case of interest occurs when the laser field intensity is very weak 
when compared to the spontaneous emission rate of a single atom {\Gi\ ^ 1). From 
f l53|l -f l55|) we find to second order in Gi that the steady-state density matrix in the 
triplet-singlet basis B f lTOj) takes the form 

[pf2^]B 
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2G? 



V 







3Fi2Gj 
V2Gi 



V2Gi 

1+5^12 
252 



1 



2- 

2(52 





2G2 



2G2 



; V2Gi 



: V2Gi 
2G2 



.(62) 



iFi2)" i-f^l2 (1-1^12) / 

Notice that Gf must be sufficiently small in order for fl62|) to make sense, since it can 
be seen that the population of the state |0,0) (component 4,4 above) diverges as the 
distance between the atoms tends to zero. In fact, it must occur that 



^ - 4 



1 (l-|F,,)^(l + fF,,)^ 



(63) 



in order that (1, 0|pfj|l, 0) + (0, 0|pf2^|0, 0) < 1. If one uses the bound in ([M]) it follows 
from ([62]) that 



[l,0\pt^\l,0) = - 



'12 



(0,0|pfj|0,0) = ^ + 



9-^12 



2 l+ilFuY 



— > as |ri — — !■ 0" 
— > 1 as Iri — — )■ 0"* 



(64) 



Hence, for a weak driving field and small distance between the atoms, the level |0, 0) is 
much more populated than the level |1,0). It must be kept in mind that the effects of 
the coherent dipole-dipole interaction hdd could modify deeply the results of this limiting 
case. 






0.5 1.0 1.5 2.0 2.571 



(b) 



Gi 

1.95 2.00 2.05 2.10 yi 

(c) 

Figure 1: (Color online) Case G2 = 0: (|Tal) and (flbl) show the steady-state concurrence 
C^'^ of the two atoms as a function of Gi = Gi/71 when = |doip and d\ = 0, 
respectively. Results are shown when the distance between the two atoms is A^i/lOO (red- 
sohd line), Aa/4 (blue-dashed line), and (black-dotted line) with \a the wavelength 
associated with the atomic transition. ( ITcl) shows a close-up of the entanglement sudden 
death when the distance between the two atoms is Aa/4. This is shown for the cases 
d^ = |doip (blue-sohd line) and d\ = (red-dashed line). 



We now discuss the degree of entanglement of pfg as measured by the concurrence 
C^^, figure ([I]). We observe that as the intensity of the laser field increases (Gi 
increases , increases, takes on a maximum value, then decreases and finally dies 
abruptly, figure ([T]). This sudden death of entanglement as a function of the laser 
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field intensity is found for most values of the distance between the atoms. Notice that 
C^'^ takes on large values only when the distance between the atoms is less than the 
wavelength associated with the atomic transition and when the intensity of the laser 
field is such that Gi = 6*1/71 < 1. This reflects the fact that the atoms can become 
entangled by exchanging spontaneously emitted photons only when they are not so 
far apart. Notice that in the region where the atoms are very close the dipole-dipole 
interaction /i^^ could be important. 

For the model under consideration, to understand the behaviour of the steady-state 
concurrence we consider the case of a weak driving field. The concurrence of the density 
operator in (l62ll is easily calculated to be 

C"""^ = ^|Fi2|-( (0,0|pfj|0,0) + (l,0|pfj|l,0)) , (65) 

to second order in Gi. Hence, we observe that for weak driving fields the steady- 
state concurrence is determined by the populations of the states |0,0) and |1,0). If, in 
addition, the atoms are close, from ( IMl) we observe that the population of the level |0, 0) 
dominates over that of the state |1, 0). Hence, the population of |0, 0) and the function 
F12 determine the behaviour of the steady-state concurrence. This can be observed in 
figure ([2]) where it is seen that the population of |0,0) is responsible to a large extent 
of the behaviour of the concurrence not only for weak laser field intensities. From that 
figure we observe that, when the distance between the atoms is less than a wavelength 
associated with the atomic transition, the concurrence grows with the population of the 
level |0,0) and is maximized when the latter is maximized. As the population of the 
state |0, 0) decreases to its asymptotic value (that is, when Gi — >■ -l-oo), the concurrence 
decreases to zero. Notice that this behaviour is not exhibited when the two atoms 
are separated by a wavelength of the atomic transition, figure fl2c|) . The reason for 
this can be explained using and fl65|) for a weak driving field. When the atoms 
are more separated, F12 tends to and we find from ( l64l) that the populations of the 
levels |0,0) and |1,0) tend to be the same. Therefore, from (!65|) we observe that the 
concurrence is determined by both of these populations and not only one of them as 
in the case when the atoms are close. Why the concurrence disappears when the laser 
field becomes very intense can be explained by the asymptotic X-form of the density 
matrix in the product state basis, (l59l) . There we observe that the laser field has reached 
such a high magnitude that many coherences have decreased to zero and that the atomic 
transitions are saturated. The latter is taken to mean that the populations have reached 
their asymptotic value. Furthermore, the concurrence of this particular X-state can be 
easily calculated to be zero. 

Also note that the decay of C'^^ with the distance between the two atoms is slower 
when d\ = (that is, the dipole moments of the two atoms are parallel to the line 
that joins them), and that C^^ exhibits an oscillatory behaviour as a function of this 
distance that is more noticeable when = |doip, figures ([1]) and ([3]). This behaviour 
is easily explained in the case of a weak laser field intensity {Gi is small). From ( 165|1 we 
see that C^^ depends on IF12I which is an oscillatory function of the distance between 
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the two atoms, fl48l) . It can be seen that the F12 has oscillations of greater amplitude 
when d\ = jdoip and that these oscillations are smoothed out when = 0. This is 
the reason why C^'^ has a greater oscillatory behaviour when d\ = |doip than in the 
case where (i^ =0. 

c' 
0.5 

0.4 

0.3 

0.2 

0.1 

0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 0.05 0.10 0.15 0.20 0.25 

(a) (b) (c) 

Figure 2: (Color online) Case G2 = 0: (I2al) and pbl) show the steady-state concurrence 
C^^ of the two atoms as a function of the population poo of the antisymmetric state 
|0,0) when d^ = |doip and d^ = 0, respectively. Results are shown when the distance 
between the two atoms is A^/100 (red-sohd line), Aa/4 (blue-dashed line), and 
(black-dotted line) with Aa the wavelength associated with the atomic transition. ( I2c|l 
shows a close-up when the distance between the two atoms is Aa- This is shown for the 
cases d\ = |doip (black-solid line) and d^ = (black-dashed line). 






02468 10 02468 10 



(a) (b) 

Figure 3: (Color online) Case G2 = 0: (!3al) and ( l3b1) show contour plots of the steady- 
state concurrence C^'^ as a function of x = uaI^i — r2\/c (horizontal axis) and Gi 
(vertical axis) for the cases d^ = |doip and d\ = 0, respectively. Contours for 
G^^ = 0.1, 0.05, 0.03, 0.01 are drawn. 

In the preceding paragraphs we have established that for any distance between the 
two atoms the entanglement present in the system disappears as soon as the laser field 
intensity is high enough (that is, Gi has achieved a sufficiently high value). So it is 
natural to ask whether all quantum correlations disappear also. Since pf^ has the form 
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of an X-state in the limit Gi +00 fl59p . we can calculate numerically very easily with 
f l26|l -( l28|) the quantum discord of pfj when measurements are carried out on one of the 
two atoms. 

We will first consider the case in which the measurements are made on the atom 
at position r2. It is found that the steady-state quantum mutual information /(pf^) is 
given by: 



-4+ y/ y 



log2(-4 + y^) 



9 P2 
2-^12 



l0g2 ( ^i^l'2 



16 + 

y 



log2 ( 16 + "-Fl^ 



9 

2^ 



where ?/ = 16 + OFfg- calculate numerically the right quantum discord. 
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Figure 4: (Color online) Case G2 = 0: (jla|l and fjlbl) show the steady-state quantum 



mutual information /"^^ (red-solid line), classical correlations C2^^'^^ (blue-dashed lined), 
and quantum discord D2 (black-dot-dashed line) as a function of x = w^|ri — r2|/c 
in the limit — )■ +00. Here measurements are performed on the atom at position r2. 
lall illustrates the case cP, = IdoiP, while (l4b1) shows the case ci^ = 0. 



Figure (jl]) shows the quantum mutual information, right classical correlations, and 
right quantum discord as a function of x = w^lri — r2|/c. The first observation is that 
pfj has both quantum and classical correlations even though the entanglement is long 
gone (from figure ([1]) we see that it died abruptly around Gi = 2). In particular, the 
quantum correlations are larger than the classical ones except for x < 2. Moreover, 
notice that the three correlations oscillate as a function of the distance between the 
atoms and that this oscillatory behaviour is much more pronounced when = |doip. 
Also, the quantum discord is exactly equal to zero for certain distances. The oscillatory 
behaviour and the exact location of the discrete zeros of the quantum discord -D^(pf^) 
are determined by F12 (see below for a transparent justification of this). Finally, notice 
that D^ipu) = HPi!) = Cf{pf2) = when |ri - r2| -> +00. Hence, pf^ has no 
correlations at all when the atoms are far apart. This is expected, since the atoms do 
not interact with each other through the reservoir in this case. 

It is of great interest to compare the right quantum discord with the geometric 
measure of quantum discord given in (l32l) to see how well the latter estimates the exact 
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result. Using that equation, the right geometric discord D^\pf2) is easily calculated 
to be 



;i6 + 9^2^)2 



(67) 
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Figure 5: (Color online) Case G2 = 0: ( l5a|) and ( l5b|) show the steady-state exact 
quantum discord -D^(pf;J^) (blue-dashed line) and the geometric measure of quantum 
discord D^\pf2) (red-sohd line) as a function of x = cjaIti — r2|/c in the limit 
Gi — )■ +00. Measurements are performed on the atom at position r2. fISal) illustrates 
the case = |doip, while (ISbll shows the case d\ = 0. 

From figure ([5]) we observe that there is good qualitative agreement between the 
estimate D^\pf2) and the exact -D^(pf^). In particular, they have the same zeros. 
Moreover, we note from fl^7|) that both the oscillatory behaviour and the zeros of 
^^\pi2) S'i'e due to F12. Hence, it follows that the discrete zeros of the exact quantum 
discord -D^(pfl') are exactly those of F12, (H5]l . Also note that the decrease in the 
amplitude of the oscillations of -D^(Pi^) when = is due to the fact that the 
oscillations of F12 are greatly decreased when d\ = 0. Nevertheless, notice that the 
exact quantum discord is considerably larger than that given in (167|) . so, in general, 
the quantitative agreement between /^^(pf^J") and Df\pf2) is not so good (except in 
determining the zeros of -D^(pf^)). 

The geometric discord in fl32l) allows us to calculate easily an estimate of the 
quantum discord for the general case flS^ - flSS]) . Figure (jH]) shows that for a fixed distance 
between the atoms the geometric discord quickly achieves the stationary value given in 
(!67|) . Notice that the decay of the geometric discord with increasing distance between 
the atoms is slower when li^ = 0. This is seen more easily in the contour plots of the 
geometric discord shown in figure ([7]). In these figures it is also observed that, as a 
function of the distance between the atoms, the geometric discord has oscillations of a 
greater magnitude when d\ = |doip. We saw above the these are due to F12. 

We now consider the case in which the measurements are performed on the atom 
at position ri. In the limiting case where IGil — ?■ +00, pf^ takes the form of an X- 
state fl59|) . The quantum mutual information /(pf^) is exactly the same as before (this 
quantity is independent of who is being measured), but we find (numerically) that the 
quantum discord D^{pi2) is now zero and the classical correlations (pf^) are equal 
to /(pfj). Hence, all information of the atom at position T2 contained in the correlations 
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Figure 6: (Color online) Case G2 = 0: f l6a|) and f l6b|) show the steady-state geometric 
discord D^^^'^ as a function of Gi = 6*1/71 for a distance between the two atoms equal 
to A^/100 (red-solid hne), Aa/4 (blue-dashed line), Xa (black-dot-dashed line) where 
Xa is the wavelength associated with the atomic transition. (J6al) illustrates the case 



doil , while (l6b|) shows the case d]^ = 0. 
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Figure 7: (Color online) Case G2 = 0: ( 17al) and ( [7b|) show contour plots of the steady- 
state geometric discord Df\pf2) as a function of x = waIi"! — ^■2|/c (horizontal axis) 
and Gi (vertical axis) for the cases d\ = |doip and d^ = 0, respectively. Contours for 
D2\pfl) = 0.035,0.025,0.015,0.005 are drawn. 



of pfl" can be extracted by measuring the state of the atom at ri [6]. This contrasts 
greatly with the case considered above, that is, with the case were the measurements 
are performed on the atom at r2. There we found that pfj presents non-negligible 
quantum discord. Thus, when measurements are performed only on the atom at r2, pfj" 
is perturbed and there is information that cannot be extracted [6]. Furthermore, we 
have here another example were both the quantum discord and the classical correlations 
are not symmetric: D^{pf:[) D'^{pf^) and Gf{pf^) ^ G^^pf^). 

Again the use of the geometric discord fl32l) allows us to give an estimate of the 
quantum discord in the general case (I53l)-(l55l). From figure ([H]) we note that the 
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geometric discord as a function of Gi is smooth by parts, since it exhibits edges. Also 
notice that it tends to zero as Gi increases, a result in accordance with that obtained 
above Df{pfl) = in the limit Gi ^ 0. Contour plots of Df\pf'^) are also shown in 
figure dH). Notice that the L'f^(pfJ) decays more slowly function of the distance 
between atoms for d\ = {] when compared to d\ = |doip, and that it is drastically 
different to D^\pf2) (compare figures ([7]) and Q). 
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Figure 8: (Color online) Case G2 = 0: (IHal) and (l8bl) illustrate the steady-state geometric 
discord D^'^^'^ as a function of Gi = 6*1/71 when the distance between the atoms 
is Xa/100 (red-solid line), Ayi/4 (blue-dashed line), and (black-dot-dashed line) 
where is the wavelength associated with the atomic transition. f l8a|) shows the case 
= |doip, while (poj) illustrates the case d']_ = 0. 
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Figure 9: (Color online) Case G2 = 0: ( l9all and ( l9b1) show contour plots of the 
geometric discord Di'\pi2) as a function of x = cjaIti — i"2|/c (horizontal axis) and 
of Gi = Gi/71 (vertical axis) for the cases d^ = |doip and d^ = 0, respectively. 
Contours for D'y \pfl) = 0.05,0.03,0.01 are drawn. 



We now turn to discuss the degree of mixed-ness of the pfj, figure (fTOj) . We observe 
that as the intensity of the laser field is increased (Gi increases), increases until it 
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(6^ 



acquires an asymptotic value which is easily calculated from fl59p : 

qST ^3 A 

^ 4 16 + 9F22 ■ 

Notice that the system is never in a maximum mixed state (as measured by the linear 
entropy) and that the maximum value 11/20 of Sf^ occurs when the atoms are very 
close. Also, the linear entropy takes an asymptotic value of 1/2 when the atoms are far 
apart. This behaviour is understandable since the dissipative interaction between the 
two atoms by means of the reservoir is strongest when the atoms are close together, and 
only one atom is affected by the laser field when the distance between the two atoms is 
large. Finally, we note that the variations of S'f^ are smoothed out when (f[_ = 0, figure 
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Figure 10: (Color online) Case G2 = 0: fllOal) and (llObI) illustrate the steady-state linear 
entropy S^^ as a function of Gi = 6*1/71 for = |(ioiP and d^ = 0, respectively. 
Results are shown when the distance between the two atoms is A^^/lOO (red-solid 
line), Aa/4 (blue-dashed line), and Xa (black-dotted line) where Xa is the wavelength 
associated with the atomic transition. 



To end this Section we discuss the relation between the entanglement between 
the two atoms as measured by the concurrence C^"^ and the degree of mixed-ness as 
measured by the linear entropy S^'^. We observe that as state becomes more mixed, the 
degree of entanglement increases until Sf^ reaches a certain value, after which the degree 
of entanglement decays to zero, figure filial) . This behaviour corresponds to the fact that 
as the intensity of the laser field increases, the state of the two atoms becomes entangled 
and mixed due to the exchange of spontaneously emitted photons. Nevertheless, once 
the intensity of the electric field is strong enough to drive the populations of the triplet- 
singlet basis ( [T0|) near their asymptotic values, the entanglement begins to decrease 
while the state of the system still becomes more mixed. We also note that the behaviour 
shown in figure filial) also illustrates the known facts that a mixed state of two qubits 
cannot contain an arbitrary amount of entanglement, and that the more mixed the state 
becomes, the less entanglement it has |l3j. 
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Figure 11: (Color online) The figures (11 lap and (lllbl) illustrate the steady-state 
concurrence C^^ as a function of the steady-state linear entropy Sf^ for the cases 
G2 = and Gi = G2, respectively. In these figures d\ = |doip. 



5. G(ri) = G{r2) 



In this section we consider that the atoms are fixed at positions where the classical 
electric field fj^Tl) has the same value. Hence, Gi = G2 which can be varied by changing 
the intensity of the laser field (see the definition of G{r) following fHT]) ). The value of the 
distance |ri — between the two atoms can also be varied independently. For example, 
such case would occur in the standing wave configuration with g{rj) = cos(kj^ ■ Vj) if 
the two atoms are located on the same line perpendicular to the wave vector k/, of the 
standing wave. 

The solution pf J of (l52|) has the following matrix representation in the triplet- singlet 
basis B ([101) : 
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^^Lt^Gi (1 + |Fi2 + AGfi 
AG\ + 2Gl 



i(l + f^i2) 







K2 = lQG\ + AGl + 









(69) 



AG\ I 



l + -F,2 



and Gi = Gi/71 again. We note that this solution is valid only if the distance 
between the two atoms is not zero. If the |ri — =0, then the solution is 

pf^ = (l-Poo)pf2^i + Poo|0,0)(0,0| , 

where Pqo is the population of the state |0,0) in the initial state, and pf2i is 



(70) 
|ri - r2| 

(71) 
obtained 
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from (1^ by replacing ^2 and (0, 0|pfj|0, 0) by 

Kg = UGt + 46-2 + 1 , (72) 

and 0, respectively, and by taking the limit |ri — — 0^ in the resulting expression. 

Note that some coherences depend explicitly on time; hence, in the IP we are 
working, a steady-state solution exists only for the resonant case = 0). The steady- 
state solution in (1691) and (I7T1) would also be approximately valid for small enough 
detuning, for instance I^l^I <^ 1. In the following we will restrict to the case where the 
steady-state solution exists, that is, we will take 5l = 0. Moreover, we will consider only 
the case |ri — 7^ 0, although we will make some remarks about the case |ri — r2| =0. 

In complete analogy to the case G2 = 0, it can be shown that, given any initial 
state pi2(0), the density operator pi2(t) of the two atoms tends to the state (l69l) above 
if the distance between the two atoms is not zero and 5l = 0. 

For the case |ri — r2| =0 and (5l = the argument is a bit different because the 
zero eigenvalue has both algebraic and geometric multiplicity equal to two, and the rest 
of the eigenvalues have negative real part. In a similar notation as before it follows that 

x(t) — i- ciVi + C2V2 if t — )• +00. (73) 

Now it can be shown that eig (the vector which has a 1 in the 16-th component and 
in the rest) and the vector associated with pfj^^ in ( ITTi) are elements of the null space 
of A, and that eig is orthogonal to the rest of the eigenvectors of A. Hence, vi and V2 
can be chosen to be those vectors, and ci and £2 must be respectively equal to Pqo and 
1 — Pqo, where Pqo is the population of |0, 0) in the initial state pi2(0). Thus any initial 
state Pi2(0) tends to the density operator in flTTl) when |ri — r2| = and 5l = 0. 

We immediately observe from ( 169|) that the effect of the collective interaction with 
the reservoir is much smaller than in the case G2 = studied above. This is due to the 
fact that the term F12 appears in combinations in which it rapidly becomes negligible 
when Gi increases (see (l69l) and flTOj) and the matrix elements there). Also notice that 
the collective interaction with the reservoir has an opposite effect in the populations of 
the triplet-singlet basis B ( ITOi) when compared with the case where G2 = 0. In the latter 
the effect is to increase the populations of the states |1,1) and |0,0) and to decrease 
the populations of the states |1,0) and |1, —1) when the atoms are close. This is easily 
seen in the high laser field intensity limit (I6OI) . On the other hand, in the case Gi — G2 
the effect is to decrease all the populations of the states of the triplet-singlet basis B 
( fTOj) . except for the state |1, —1). Nevertheless, the effect is small, since the terms in Gi 
quickly overwhelm those with F12 as the laser field intensity increases (Gi = G2 grows). 

There are several limiting cases of interest. First, if Gi — )■ 0, then it is seen from 
fl69|) that pfj —7- 1|. This result is expected since the atoms are located at 

different positions and, without the driving field, the reservoir ultimately leaves the two 
atoms in their respective ground states. Notice that this result would have not been 
obtained if the atoms where exactly in the same position |ri — r2| =0 because in that 
case the antisymmetric state |0,0) is invulnerable to dissipation. 
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When the laser field intensity is weak (|Gi| = 16*21 ^ 1) we find from fl69p that the 
matrix representation of pf^ in the triplet-singlet basis B (11 01) takes the form 



V 
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2+3Fi2 
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0/ 



(74) 



to second order in G\. Notice that G\ must be sufficiently small in order for (17^ to be 
valid. In fact, it must occur that 



G\ < 



(2 + 3F] 



12 



32 



(75) 



to have a population of the state |1,0) less than or equal to 1. Hence, the symmetric 
state |1, 0) is much more populated than the antisymmetric state |0, 0) for a weak laser 
field intensity {Gi = G2 is small). Notice that this contrasts greatly with the case 
G2 = (see (l64l)). From (17^ it is also observed that the collective interaction with the 
reservoir tends to decrease the population of the state |1, 0) and to increase that of the 
state |1, —1) when the atoms are close. It is important to keep in mind that for very 
small distances between the atoms the effects of the coherent dipole-dipole interaction 
not incorporated in the model under consideration can be very important. 

Another limiting case of more interest occurs when the electric field is very intense 
(\Gi\ = IGal +00). From ([69]) we note that pf:[ = (1/4)1 with I the identity 
operator. This behaviour is expected, since for high field intensities the atom-laser field 
interaction overwhelms the atom-atom interaction through the reservoir and each atom 
approximately interacts independently with the laser field. Also notice that in the limit 
of high field intensity, the populations of the states |1, 1) and |0,0) take their absolute 
maximum value of 1/4, while that of the state |1, — 1) takes its absolute minimum value 
of 1/4. 
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Figure 12: (Color online) Case Gi = G2'- (I12al) and (112b|) show the steady-state 
concurrence C^'^ as a function of Gi = 6*1/71 for the cases = |doip and d\ = 0, 
respectively. Results are shown for distances between the two atoms equal to Ayi/100 
(red-solid line), Xa/4: (blue-dashed line), and Xa (black-dotted line) where Xa is the 
wavelength associated with the atomic transition. Figure f ll2cp shows a close-up of the 
cases d^ = |doip (black-sohd line) and d\ = (black-dashed line) when the distance 
between atoms is A^^. 
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We now discuss the degree of entanglement of pf^ as measured by the concurrence 
C^^, fi gure f|T2|) . We observe that, as the intensity of the laser field increases (Gi = G2 
increases), C^'^ increases, takes on a maximum value, then decreases and finally dies 
abruptly, figures (fT2|) and (fl^ . This sudden death of entanglement as a function of the 
laser field intensity is found for most values of the distance between the atoms. Also 
notice that C^'^ takes its largest values when the distance between the atoms is less 
than a wavelength associated with the atomic transition and when the intensity of the 
laser field is small, although in this case the largest values of C^'^ are much smaller 
than those of the case G2 = 0, figures ([1]) and ( !T2|) . It thus appears that significant 
entanglement between the two atoms cannot be generated when they interact only 
through the reservoir and if Gi = G2 (recall that we are only considering the dissipative 
interaction). This is physically understandable when the the laser field intensity is large, 
since each atom interacts independently (approximately) with the laser field due to the 
fact that this interaction overwhelms the atom-atom interaction through the reservoir. 
This is confirmed by the weak dependence of pfj on F12 in (l69l) and ( ITTll . On the other 
extreme, when the laser field intensity is small, from ( 17^ we find to second order in 
Gi = G2 that the concurrence of pfj takes the form 

max<0, -^Fi2(l,0|pf2'^|l,0)l . (76) 



C 



ST 



Now F12 is an oscillatory function bounded above by 2/3 and bounded below by —0.2237 
which tends to as the distance between atoms increases. Hence, the concurrence to 
second order in Gi in (1761) will be non-zero only for certain values of the distance between 
the two atoms. The equation in fl76|) might lead one to think that the behaviour of 
the concurrence is determined to a great extent by the population of the state |1,0). 
From figure (fT3|) we see that this is not the case, since the population still increases 
considerably after the concurrence has already vanished. 
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Figure 13: (Color online) Case Gi = G2: (I13al) and (]13b|) show the steady-state 
concurrence G^'^ as a function of the population pio of the state |1,0) for the cases 
d^ = |doip and d^ = 0, respectively. Results are shown for distances between the two 
atoms equal to Aa/100 (red-solid hne), (blue-dashed line), and (black-dotted 

line) where A^ is the wavelength associated with the atomic transition. Figure f ll3cp 



shows a close-up of the cases d"] 



101 1 



(black-solid line) and d^ = (black-dashed 



line) when the distance between atoms is A^- 



Also observe that G'^^ decreases its value when the distance between the atoms 
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increases, and that this decay with distance is slower when d\ = (that is, the dipoles 
of the two atoms are parallel to the line joining them), figures f|T2|) and f|T^ . Finally 
notice that C^^ exhibits an oscillatory behaviour when = |doip (that is, the dipoles 
of the two atoms are orthogonal to the line joining them), and that this phenomenon 
is smoothed out when d^ = 0. Again, the reason for this is the function F12 that has 
oscillations of greater amplitude when d\ = |doip. 
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Figure 14: (Color online) Case Gi = G2: (I14ap and (114bp show contour plots of the 
steady-state concurrence C^^ as a function of x = uaI^i — r2\/c (horizontal axis) and 
Gi (vertical axis) for the cases d^ = |doip and d^ = 0, respectively. Contours for 
= 0.07, 0.05, 0.03, 0.01 are drawn. 



As before, we will now study the behaviour of the quantum discord of 
function of both the distance between the atoms and of Gi. We first observe that in 
this case -Of (pf^) = -D^(pf;f ), since the two atoms are in equivalent positions, that is, 
both atoms are in a position where Gi = (j2. Note that this contrasts with the case 
considered in the previous section where the atoms were in non-equivalent positions and, 
as a result, D?{pf^) ^ D^{pfl). 

We observed above that pf^ = (1/4)1 in the limit Gi — )■ +00. Hence, we 
conclude from ([1]), ([7]), and ([8]) that pf^ has no correlations in this limit, that is, 
= I{pfl) = Cfipfl) = Dfipf^) = (j = 1,2). This is explained intuitively 
by the fact that in the high laser field intensity limit (that is, \Gi\ — )■ +00) each 
atom interacts independently (approximately) with the laser field and the atoms are 
in different positions (if the atoms were exactly at the same position, the quantum 
discord would have to be calculated from flTTl) ). 

Using ( 132|) we can study the geometric discord Df\pf2) in the general case ( l69l) . 
From figure (fTSll we observe that the D^^\pf2) tends to zero as Gi increases, a result 
expected due to the fact that (pf^) = in the limit \Gi\ — )■ +00. Moreover, notice 
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that D^i\pi2) decreases more slowly as the distance between atoms increases when 
d\ = 0. This is clearly seen in figure ( |T6l) where it is also observed that Df\pf2) 
has oscillations of a greater magnitude when d\ = |doip. As before, the oscillatory 
behaviour of -Dp\pi2^) is due to F12, a result which is known due to the analytic but 
unwieldy expression obtained from the use of fl32l) . 

From the results presented above it thus appears that driving the two atoms by the 
laser field is ineffective in creating quantum correlations, since both the entanglement 
and the quantum discord are considerably smaller when compared to the case where 
only one atom is being driven (the case G2 = considered above). This decrease of 
correlations is even more dramatic in the limit of high laser field intensity (IGil — ?■ +00), 
because it was shown that all correlations disappear in the case Gi — G2- 
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Figure 15: (Color online) Case Gi = G2: (11 Sal) and (115b|) show the steady-state 
geometric discord D^^^^ as a function of Gi = 6*1/71 for the cases d^ = |doip and 
d^ = 0, respectively. Results are shown for distances between the two atoms equal to 
Aa/100 (red-solid line), (blue-dashed line), and Xa (black-dotted line), where Xa 

is the wavelength associated with the atomic transition. 



We now turn to discuss the degree of mixed-ness of pf^J", figure (fT7|) . As the laser 
field intensity increases (Gi grows), Sf^ increases until it acquires an asymptotic value 
of 3/4. In fact, we observed above that pf^ = (1/4)1 when IGil — t- +00. Hence, pf^ 
tends to a maximum mixed state when the laser field intensity increases. Again, this 
behaviour is expected, since for high field intensities each atom approximately interacts 
independently with the laser field. Finally, notice that Sf^ depends weakly on the 
distance between the two atoms and on the value of d\, figure (fTTj) . This is expected, 
since it was observed above that pf^ depends weakly on F] 



12- 



To end this Section we discuss the relation of the steady-state entanglement between 
the two atoms as measured by the concurrence G^^ with the steady-state degree of 
mixed-ness as measured by the hnear entropy Sf^. We observe that, as Sf^ becomes 
larger, the degree of entanglement increases until a maximum value. After this the 
degree of entanglement decays to zero, figure flllb[) . This behaviour corresponds to 
the fact that, as the intensity of the electric field increases, the state of the two atoms 
becomes entangled and mixed due to the exchange of spontaneously emitted photons. 
Nevertheless, once the intensity of the electric field is strong enough to dominate over 
the interaction between the atoms by means of the reservoir, the entanglement begins 
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Figure 16: (Color online) Case Gi = G2- (llGap shows a contour plot of the geometric 
discord D^^^^ of function of a; = UA\ri — r2\/c (horizontal axis) and Gi (vertical 

axis) for the case ci^. = |doiP- Contours for D^^^"^ = 0.005,0.003,0.001 are drawn. 
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Figure 17: (Color online) Case Gi = G2: fll7al) and fll7bl) show the steady-state linear 
entropy S^'^ as a function of Gi = 6*1/71 for the cases d\ = |doip and d\ = 0, 
respectively. Results are shown for distances between the two atoms equal to \a/100 
(red-solid line), (blue-dashed line), and \a (black-dotted line) where Xa is the 

wavelength associated with the atomic transition. 

to decrease while the state of the system still becomes more mixed. Figure (lllbp 
also illustrates the known facts that a mixed state of two qubits cannot contain an 
arbitrary amount of entanglement, and that the more mixed the state becomes, the less 
entanglement it has [13j . 

6. Conclusions 

In this article we considered two two-level atoms (qubits) fixed at different positions, 
driven by a monochromatic laser field, and interacting collectively with the vacuum 
electromagnetic field. A Born-Markov-secular master equation was used to describe the 
dynamics of the two atoms and their steady-state was studied for two configurations. 
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In the first, one atom was located in a region where the driving electric field is zero 
while the other was in a position where it is not zero. In the second configuration both 
atoms were located at different positions where the driving electric field was the same. 
We refer to the first case as 6*2 = and to the second case as Gi = G2- 

We neglected the coherent dipole-dipole interaction between the two atoms that 
results from the collective interaction with the vacuum reservoir, and we only kept 
the dissipative atom-atom interaction through the reservoir. This led us to a model 
of interest in itself: two driven qubits at different fixed positions interacting with a 
vacuum reservoir whose Born-Markov-secular master equation is (150|) . Furthermore, 
this model allowed us to identify and understand which effects are due solely to the 
dissipative interaction between the two atoms. For this model the steady-state density 
operator exists in an interaction picture only when the laser field frequency is resonant 
with the atomic transition frequency (the resonance condition). The analytic steady- 
state density operator pf^ was obtained for both configurations of the atoms. This 
allowed us to study several limiting cases that resulted in a much deeper understanding 
of the steady-state of the system and the mechanisms responsible for the formation of 
correlations between the two atoms. 

It was shown for both configurations that any initial state of the two atoms tends to 
the respective pfj under the resonance condition. These steady-state density operators 
are entangled and have non-zero left and right quantum discord if the laser field intensity 
is not very high. High laser field intensities turn pfj into a maximum mixed non- 
correlated state (1/4)1 in the case Gi = G2, and into a separable X-state that has 
non-zero right quantum discord (measurements are made on the un-driven atom) in the 
case G2 = 0. 

It was found that driving both atoms with the laser field (case Gi = G2) is not 
efficient for the generation of steady-state quantum correlations (entanglement and 
quantum discord), since all steady-state quantum correlations decreased their values in 
the case Gi = G2 and only the steady-state degree of mixed-ness was increased. In both 
configurations, steady-state quantum correlations (entanglement and quantum discord) 
were maximized for weak laser field intensities. Increasing the laser field intensity led 
in all cases to steady-state entanglement sudden death, but the quantum discord still 
survived and, in some cases, did not decrease to zero. We also found that the distance 
between atoms is fundamental for the build up of steady-state quantum correlations. 
The farther the atoms are apart, the less quantum correlations they can have. This is 
physically expected since the atoms can only interact by interchanging spontaneously 
emitted photons and this interaction is less effective the farther the atoms are apart. 

In general, all correlations exhibited an oscillatory behaviour that tended to zero as 
the distance between the two atoms was increased. The oscillations were more marked 
when the dipoles of the atoms were orthogonal to the line joining them and the decay 
with distance was slower when the dipoles were parallel to the aforementioned line. It 
was found that this behaviour was due to the function F12 that comes about from the 
collective interaction of the atoms with the reservoir. 
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The configurations considered also allowed us to study the behaviour of the 
quantum discord when the atoms are in non-equivalent and equivalent positions. In 
particular, a method to calculate numerically the left and right quantum discords was 
presented, and it was used in the case of a high laser field intensity. This allowed 
a direct comparison with the geometric measure of quantum discord. It was found 
that the geometric discord reproduced the behaviour of the quantum discord and there 
was good qualitative agreement between both. The geometric measure was also used 
to study the cases where the laser field intensity was not high and showed that the 
quantum discord survives long after the steady-state entanglement has disappeared. 

Due to the non-equivalent position of the atoms in the G2 = configuration, it 
was found that the right geometric discord D^2\pi2) (obtained when measurements are 
made on the un-driven atom) is very different from the left geometric discord -Di^''(PiJ) 
(obtained when measurements are made on the driven atom). D^\pf2) rapidly acquires 
a stationary non-negligible value as the intensity of the laser field is increased and has 
an oscillatory behaviour as a function of the distance between atoms with discrete 
zeros determined by the function F12. On the other hand, -^^^^(pfj) tends to zero as 
the laser field intensity is increased. Hence, the information that cannot be extracted 
from measurements only on the driven/un-driven atom is very different in the ^2 = 
configuration. This is not the case of the Gi — G2 configuration where the atoms are in 
equivalent positions. 

Analytic expressions for the steady-state populations of the eigenstates of the free 
Hamiltonian of the two atoms where derived and we established their simple behaviour 
in the limits of high and low laser field intensities. We determined that the populations 
of the antisymmetric |0, 0) and symmetric |1,0) states explained the behaviour of the 
concurrence to a great extent in both configurations. This was much more dramatic 
in the case G2 — ^ were we showed that the concurrence increased/decreased with the 
population of the antisymmetric state almost linearly when the distance between the 
two atoms was less than a wavelength of the atomic transition. 

In this work we have characterized the quantum correlations between the two 
atoms when they interact only through the dissipative collective interaction. It is very 
interesting to study the effects of the dipole-dipole interaction between the two atoms, 
not taken into consideration in this work, to see what new dynamics will be introduced 
or if the system will behave similarly. This will be investigated in future work. 
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